Collective Modes in Two-band Superconductors. 
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We analyze collective modes in two-band superconductors in the dirty limit. It is shown that these 
modes exist at all temperatures T below T c provided the frequency of the modes is higher than the 
inelastic scattering rate and lower than the energy gaps A a ^- At low temperatures these modes 
are related to counterphase oscillations of the condensate currents in each band. The spectrum of 
the collective oscillations is similar to the spectrum of the Josephson "plasma" modes in a tunnel 
Josephson junction but the velocity of the mode propagation in the case under consideration is 
much lower. At higher temperatures (Aj, < T < T c ) the spectrum consists of two branches. One of 
them is gapless (sound-like) and the second one has a threshhold that depends on coupling between 
the bands. We formulate the conditions under which both types of collective modes can exist. 
The spectrum of the collective modes can be determined by measuring the I-V characteristics of a 
Josephson junction in a way as it was done by Carlson and Goldman [12 ] . 

PACS numbers: 74.20.De, 74.25. Ha, 74.25.-q 



a 

o 

O i I. INTRODUCTION 

U 

The conventional BCS theory has been developed for a single-band metal with an attractive interaction between 
electrons with opposite spins and momenta (s-wave, singlet pairing) | lj and describes well most low-T c superconductors. 
The universality of the description of the superconductors is a consequence of an assumption about a simple shape of 
the Fermi surface. 

However, some superconducting materials have a rather complicated band structure. For example, there is a 
consensus that the recently discovered new superconductor MgB^ (T c w 40 K) is a two-band superconductor [2]. In 
contrast to conventional BCS superconductors, two-band superconductors may have two different order parameters 
A a ,fc (we label the bands by subscripts a and b) and an additional degree of freedom - the phase difference between 
the order parameters: ip = Xa — Xb- This is a new variable that has to be accounted for in the proper theory of such a 
i— i , superconductivity. Naturally, one can expect new phenomena in two-band superconductors related to this new degree 
of freedom. 

One of the examples of this kind are p— phase solitons predicted in On the basis of the Ginzburg-Landau (GL) 
functional it was shown that these solitons are described by the sine-Gordon equation. A Gedanken-experiment that 
■^j- | allows the observation of these solitons was suggested in [J] . The authors considered a two-band superconductor in 
fT} . contact with a normal metal N. Using a version of the time-dependent GL equation, they analyzed a current through 
the S /N interface and showed that under certain conditions the phase solitons may be created in the superconductor. 
, Equation describing dynamics of these solitons is similar to the one for a dissipative Josephson junction. In Ref. [5| an 
analogy between the Ginzburg-Landau functional for two-band superconductors and for some models in the particle 
f^*- , physics (an extended version of Faddeev's nonlinear a model) was used and, on this basis, topologically different 
> vortices in two-gap superconductors have been predicted. In Ref.jy] the critical current in these superconductors was 
calculated. 

Another effect that can arise in two-band superconductors is weakly damped oscillations of the phase difference ip of 
the order parameters or, in other words, collective modes (CMs) related to oscillations of the phase <p in space and time. 
It is known that in single-band superconductors CMs can exist only at temperatures close to the critical temperature 
T c (see reviews @, @). In these modes the condensate current js and the quasiparticle current j q oscillate. Because 
the variation of the total current density Sj should be zero due to the quasineutrality condition (Sj = Sjs + Sj q = 0) , 
the oscillations of Sjs are accompanied by counter-phase oscillations of the quasiparticle current Sj q ~ oE. The phase 
of the order parameter also oscillates but the amplitude of the order parameter A remains constant. These phase 
modes have an acoustic spectrum uj ~ kv cm and exist in a sub-gap region (A 2 /T < lu << A, 1/t; where r is the 
elastic scattering time) in imp ure superconductors. The CMs have been observed by Carlson and Goldman [l2j and 
explained theoretically in [lfj, [llj . 

It is clear that CMs cannot exist in conventional superconductors at zero temperature because in this limit only 
one degree of freedom, namely, the condensate current js exists. Any oscillations of the current density js with not 
very high frequencies would lead to violation of the charge neutrality. In contrast, in two-band superconductors, even 
at zero temperature T there are two degrees of freedom: the condensate currents jsa,b in each band that can oscillate 
in counterphase such that the total current density remains constant: Sjs(r,t) = Sjsai^jt) + Sjsb(r,t) = 0. This is 
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similar to what happens in layered superconductors [l3T ] where the condensate currents in different layers oscillate in 
counterphase. The CMs in the two-band superconductors at zero (or low) temperature arc related to oscillations of 
the phase difference tp. These CMs are similar to the Leggett mode that can be excited in superfluid He 3 fl4l|. 

Theoretically, the Leggett- type CMs in two-band superconductors at zero temperature were studied in [la ], and 
the influence of the CMs on the Josephson effect in S1/I/S2 was analyzed in [l(| (here Si 2 is a single- and two-band 
superconductor). This Leggett-type CM has a spectrum 



n 2 = n 2 a + v 2 L k 2 , (1) 

which is typical for the Josephson tunnel junction, where Oo is a threshold frequency, vl is a velocity of this CM (see 
Sec. III). The similarity between the Josephson "plasma waves" and the Leggett-type CM in two-band superconductors 
is quite natural because the coupling between different bands in two-band superconductors looks like the Josephson 
coupling between superconductors in a Josephson tunnel junction. 

According to estimates carried out in Ref. [Til ], the energy corresponding to the threshold frequency flo in MgBi 
is higher than the smaller superconducting gap. This corresponds to a strong damping of the CMs in this two-band 
superconductor. Therefore it would be of interest to investigate under what conditions the gapless Carlson-Goldman 
CMs can propagate in two-band superconductors and this is the subject of the present paper. 

The paper is organized as follows. In Sec. II we formulate a model of two-band superconductors that will be used 
in calculations and present microscopic equations for quasiclassical Green's functions. These equations determine the 
spatial and temporal behavior of the retarded (advanced) Green's functions g^ A \t,t';r), as well as of the Keldysh 
function gf (t,t';r). The functions gf" A ^ and gf are matrices in the particle-hole (Gor'kov-Nambu) space. Using 

these equations we find a linear response 5gf^ and Sgf to small perturbations of the electric field and condensate 
velocities in both the bands. This give us possibility to find the spectrum of the CMs at arbitrary temperatures. Note 
that the spectrum of CMs cannot be obtained from the generalized GL equations used in [4j. In Sec. IV we analyze a 
method that may enable one to observe the CMs. 



II. MODEL AND BASIC EQUATIONS 



We start our discussion considering a simple model of a two-band superconductor and deriving microscopic equations 
for quasiclassical Green's functions. Using these equations we can obtain equations for macroscopic quantities and 
calculate the spectrum of CMs. We restrict ourselves with the dirty limit assuming that the elastic impurity scattering 
time r is sufficiently small: A a ,f,r << 1. In the equilibrium case these equations can be reduced to an Usadel-like 
equation used in [171 ] . 

The Hamiltonian of the considered two-band superconductor has the form (see for example [HI]) 



H = H a + H b - { V i,k(q)i>tp+ q ^t-p^k,-p>'ll)k,p'+q + CC.} , (2) 

{p,q;i,k} 

with 

Ha = ^ i^ap v a(P ~ Pa^lp + ^(p+g) [V a ,imp{q) + V a (q)]lp ap } (3) 

{p,q} 

where the third term in Eq.Q describing the electron-electron interaction leads to the superconductivity; p, q are 
momenta (strictly speaking, we must assign to p also a spin index a but we omit it here for the sake of brevity) 
and the indices {i,k} numerate the bands {a,b}. The first two terms are one-particle Hamiltonians for each band 
that include the kinetic energy s a ,b(p) = ^a,&(p — Pa,b) counted from the Fermi level and the terms describing an 
elastic impurity scattering. Beside a short-range potential V a ^ m p due to impurities, H a contains also a long-range 
self-consistent potential V a due to coulomb interaction. The Hamiltonian for the b band, H a , is obtained from H a by 
replacing subindeces a — > b. We neglect the interband impurity scattering (arguments supporting this assumption 
have been given in (l9j). 

The derivation of the equations for the quasiclassical Green's functions is carried out in a standard way [20L [2l| . 
However, the presence of the two bands makes the situation more complicated. In order to avoid unnecessary techincal 
difficulties, we make several assumptions: a) as in the BCS theory we use the mean-field approximation representing 
the product of four ip operators in the form A^-i/^ -p'Ykp'+q + p'+q' b) we neglect a change in the single- 

electron spectrum due to a possible tunneling between the bands because the strong impurity scattering assumed here 
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destroys such a change, c) we neglect terms corresponding to pairing of electrons from different bands (ipj _ p ,ipl p , +q ) 
(such a pairing was taken into account in [22j where the possibility of triplet pairing in a clean layered superconductor 
was analyzed). With these assumptions used also in previous works, we can derive a micro scop ic equation for the 
matrix quasiclassical Green's functions g in the same way as for a one-band superconductor [20L l2lj . This equation 
in the dirty limit (rA a ^ << 1) has the standard form 

- iD a V(gVg) a + i[hdg a /dt + dgjdt'h] + [A a ,g a ] + s a [A b ,g a ] - e [V{t)g - gV{t')\ = 0; , (4) 

where D a = (i>i?Z) a /3 is the diffusion coefficient in the a band, g a (r; t, t') is a 4x4 matrix depending on the coordinate 
r and two times t and t'. The elements of this matrix are the retarded (advanced) matrix Green's functions g R ( A > 
(elements (11) and (22)) and the matrix Keldysh function g K (element (12)). The parameter s a — V a b/Vb determines 
the strength of the coupling between superconducting pairing in the a and b bands (this type of pairing for the two- 
band superconductors was suggested earlier in (23|). The same equation as Eq. (T?|) is valid for the b band provided 
the subscripts are exchanged, a t» b. 

Eq.([5]) is supplemented by the normalization condition 

g a {t,ti)og a {ti,t') = 6(t-t') (5) 

and the self-consistency equation 

A a ,6 = \ a ,bf a ,b{t,t;r) (6) 

where A a .fc = {Vv) a .b and is a ,b are the coupling constant and the density of states in each band. If we wrote down 
Eq. Q for the retarded (advanced) Green's functions in the Matsubara representation, we would obtain a generalized 
Usadel equation. Such an equation was used by Koshelev and Golubov [l7j . Note that our definition of the order 
parameter differs from the one (A a )^. G used in Ref. [TtJ and the correspondence between both the definitions is given 
by {A a ) KG = A a + s a A b . 

The current density in each band is expressed in terms of the Keldysh matrix as 

j«,6(*,r) = ^/A)<j a ^Tr{T Z [9 R {tM;r)^g K {tut-r) + g^(t, i i; r)Vg A (^, t; r)] a , & } (7) 

Our aim is to find the response 5g a ,b of the system to a perturbation of the electric potential V(r, t) and ^Xa,b- 
To be more precise, we are interested in the response to perturbations of the gauge-invariant potential fi a> b and the 
condensate momentum Q a ,b 

f* a ,b = eV+ (l/2)d Xa ,b/dt] Q a , b = (l/2)[V Xa ,b - (2tt/$ )A] (8) 

where $o = hc/2e is the magnetic flux quantum. These are the responses that enter physical quantities. One more 
quantity s a [Ab,g a ] will also be considered as a perturbation with A;, and g a equal to their equilibrium values. 
In equilibrium the Keldysh function g^ h {e) equals 

= (g R (e) - g A (e)U tanh(e/3) (9) 



A 2 h 

a.b 



where g^\e) = [r 3 g a .b(e) +i f 2 f a<b (e)]^ , f$*> (e) = A^/^T = (A a , b /e)g% A \e), £j?\e) = J(I±W 
and (3 = (1/2T) . As usual, the quasiclassical functions g and / stand for the normal and condensate quasiclassical 
Green functions and t* are the Pauli matrices in Gorkov-Nambu space. The method of solution we will use is similar 
to that presented in [f| . 

First, we single out the phase \a,b using the transformation g a> b = {Ug new U^) a j,, where U a ,b = exp(if3Xa,fc/2). 
Then Eq.(f4]) for the new matrix g new acquires the form (we drop the subindex new) 

-*I>oV(dVd)a+i[^^ + ^^] + [Aa,i>a]+*o[Ai,^a]-[A*(t)d-M*')]o+- D «V ' Qaffafc, 9a]+iD a Ql[f 3 , ffa^a] = 0, 

(10) 
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After that we linearize Eq. (jTO|) with respect to perturbations Sg(t, t'\ r) = g — g eq and make the Fourier transfor- 
mations 

Sg(e,€,k) = J dtdt' exp(iet- ie't')Sg(t,t',k) (11) 

where Sg(t,t',r) ~ Sg(t, t' , k) exp(ikr). We represent the perturbations /J.(t,r) and Q(t, r) in the form: fi(t,r) ~ 
k) cxp(«kr — iflt); Q(t, r) ~ Q(fi, k) cxp(ikr — iili). 
Writing down equations for elements (11) and (22), i.e., for the matrices g R and we can obtain the expressions 
for the perturbations Sg R (e, e') of the retarded Green's functions 

e ') = TfR i 3 ^ - 9aAbC') + VatCC' - 1) - iDJsQ a (ft& - (12) 

where M* (e, e') = (ff + 0) a + ik 2 D a , (£f )„ = £*(e) and £*(e) is defined in Eq.©. 

The matrices of the perturbations Sg R (e, e') and Sg^(e, (!) arc determined by Eq. (fT2|) after the permutation of the 
subscripts: a — » 6 and i? — * A. Eq. (|12[) coincides with a corresponding equation in [8| provided the limit At << 1 is 
taken and s a is set to zero: s a = 0. 

In order to find the perturbation of the Keldysh function Sg a (e, e') = 6gff (e, e'), we represent Sg a (e, e') in the form 
of a sum of a regular ^g re& and anomalous 5g an part 

<^a(e, e') = {5g reg {e, e') + 5g an (e, e')) a (13) 

where 8g reg {t, e') — 5g R (e, e) tanh(e'/3) — tanh(e(3)Sg A (e, e ) (we drop the indeces a, b). The anomalous part is obtained 
in a way similar to that in It has the form 

(5g an (e, e')) a = tanh(£,/3) Ma tanh(£/j) {- gQ (A b - ^A^) + ^(1 - ^) - iX>„kQ (§*fa - f 3 ^)} (14) 

where M a (e, e') = (£f + 0) a + ik 2 D a . The energies e, e' in Eas. (|12H14p are equal to e = e + fi/2, e' = e - 0/2, where 
e = (e + e')/2 and f2 is the frequency of oscillations. 

Having determined the perturbations Sga ( e i e ') an d <%i(e,e')i we can readily derive equations for such macro- 
scopic quantities as fJ* a ,b{t,r),5j ai b(t,r), etc, in each band and obtain the spectrum of the CMs. 



III. MACROSCOPIC QUANTITIES. SPECTRUM OF OSCILLATIONS. 

As follows from its definition, the condensate momentum Q a ,6 obeys the equations 

dQ a ,b/dt = eE + V/i a , 6 (15) 

In order to obtain an equation for fi a ,b, we can use the self-consistency equation © written for the phases Xa,b- This 
means that terms proportional to t\ in Eq.® should be equal to zero. The variation of the current density is found 
from Eq.©. 

We consider first the case of low temperatures: T << A a b- 

a) T « A ai b. In this case the main contribution is due to the regular part Sg reg (e, e'). The anomalous part gives 
small corrections of the order fi/A^f, because we assume that £l/A a ^ << 1. 

Let us obtain an equation for fi a j, using Eq.©. Calculating the contribution from the regular part and setting 
e = iuj + Q/2 and e' = iu> — 51/2 we can transform the integration over e into a sum over the Matsubara frequencies: 

J de[5g R (e,e')t!mh{e'/3) - tanh(e/3)(5g A (e, e')] = {2m){2T) ^ 6g R (e, e) 



Substituting the perturbations 5g into Eq. © we get 
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- eosin^- u a dfj, a /dt + A a (V a D a )VQ Q = (16) 

where eg = 2(y a b/V a Vbv)A a Ab, v a = v a j(y a + v b ) is the normalized density-of-states in the a band, v = v a + Vb, and 
tp = Xa — Xb- The same equation with the interchange of indices a b is valid for the 6 band (one has to keep in 
mind that the sign in front of sin^ changes as a result of this interchange). 

Eq. (|16p is the continuity equation for the charge of Cooper pairs q a ~ v a \i a in the a band because the third term 
is proportional to the divergence of the condensate current j a . The first term in Eq. (fT6|) describes a Josephson-likc 
coupling between the bands and may be considered as a drain (ip > 0) or source (ip < 0) of Cooper pairs in the a 
band. Eq. (|16p is valid for frequencies exceeding the effective relaxation time Ti m b for a charge imbalance, Q >> l/rj m f,. 
This relaxation time is determined, in particular, by the electron-phonon and electron-electron inelastic scattering. 
If charge imbalance relaxationprocesses are taken into account, the term d/j, a /dt in Eq.(fT6|) should be replaced by 
(0/0* + 7), where 7" 1 = r lmb 0, 1 ©]. 

At low temperatures the current density j a coincides in the main approximation with the condensate current and 
equals 

ja = 7TCT Q A a Q a /e (17) 

This expression and Eqs. (fT51) , (TTB|) describe the system at low temperatures. One should add also the charge 
neutrality condition 

5j=6(j a +j b ) = 0- (18) 
Then, we can exclude all the variables except ip and obtain the equation for the phase difference 

sin ip + (d/dt + ^dip/dt - v 2 cn V 2 ip = (19) 
where S7q = 2eq/ \v a v^\ 7 = l/rj mo is a damping rate and 



v cm = [nA a A b n - A 7^ _ A ] 1/2 (20) 
D a v a A a + DbVbAb 

is velocity of the CMs or the limiting velocity of phase solitons in the two band superconductors. 

Eq. (|19|) is similar to the sine-Gordon equation for a tunnel Josephson junction but the velocity v cm is much smaller 
than the corresponding velocity of the Swihart waves in a Josephson junction. By the order of magnitude the 
velocity v cm is equal to the velocity (~ \J DA) of the Carlson-Goldman CM in ordinary superconductors. However, 
in contrast to the Carlson-Goldman CMs in single-band superconductors that are weakly damped only near T c , the 
CM described by Eq. ([T9|) exist at low temperatures. In these modes, condensate in the a band oscillates with respect 
to the condensate in the b band and these oscillations are accompanied by oscillations of the phase difference ip. The 
spectrum of the small amplitude oscillations is given by 



n 2 = n 2 Q + k 2 v 2 cm . (21) 

The threshold frequency CIq is analogous to the Josephson "plasma frequency" . As we assume that f2 is less than 
A a .fc, our consideration is valid provided V a b(v a + Vb) /2V a VbV a Vb < A /A a . This type of CMs that is analogous to 
the Leggett mode in He 3 was studied theoretically in Refs.[l5|, [l6j]- The dependence of fi(fc) is depicted in Fig. la for 
certain parameters of the model. 

Now we consider the limit of temperatures close to T c . 

b) A a ,b « T. The calculations in this case become more combersome because the anomalous Green's function 
5g a n,a{^i e')> Eq. (|14[) . also gives an essential contribution. In this case the currents j a ^b are equal to 



= a a , b [(nAi !b /2T)Q atb /e + E] (22) 

where the first term is the supercurrent and the second term is the quasiparticle current j q = aE. In the main 
approximation this current originates from the anomalous part Sg an . a (e, e'), Eq. p4p . 
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FIG. 1: a) Calculated CM spectrum for the case of low temperatures (T << A ai6 ). There is a gap flo in the CM spectrum. 
(k = f2o /vcm) b) The CM spectrum for high temperatures (T » A a ^). There are two branches of the CMs: the sound- like 
mode and the mode with a gap (analogous to the Leggett mode). The sound- like mode exists for frequencies greater than 
inverse energy relaxation time: Q > 7. 

The equation for [i a acquires the form 

l 2 sin^ + {vp) a {d/dt + 7 ) Ma - (Pp)„i£VQ a = (23) 

where Sq = (4T/7tA)£q, A = A a + A b ,p a — P A a /(A,,+A(,) and v a = ^/2D a A a is the velocity of the Carlson-Goldman 
mode in the a-band. In the limit {7, A 2 b /T} < £1 < A a , b the equation for Q a ^ is reduced in the main approximation 
to dQ a ,b/dt w V/x 0) h. The spectrum of the CMs consists of two branches determined by the roots of equation 

(fi 2 - k 2 v 2 a )(n 2 k\i) = 2iU(tf - k 2 v 2 a )/ Pb + (Q 2 - k 2 v 2 )/ Pa )} (24) 
In the long-wave limit these branches are described by the expressions 

nl - 2i 2 (l/ Pa + l/p b ) + k 2 (v 2 aPb + v 2 p a )/(p a + Pb ) (25) 

and 

n 2 = k 2 (v 2 aPa + v 2 p b )/(p a + Pb ) (26) 

Therefore one branch of the CMs has a sound- like spectrum (f^) and another one a Josephson-like spectrum (fii). 
If v a >> v b and A a >> A b (this limit corresponds to the two-band superconductor MgB 2 ), one has f2 2 = fc 2 w 2 , 
that is, the low-frequency mode coincides with the Carlson-Goldman mode in the band with a higher gap. This low 
frequency mode may be excited in such two-band superconductor as MgE>2- One can show that if HQ < A 2 /T, then 
the sound-like branch of the CMs is strongly damped. Therefore at low temperatures, the "soft" mode can hardly 
exist. Below we clarify this point in more detail. The form of the spectrum at high temperatures is shown in Fig. lb. 
Consider now the case of intermediate temperatures, 
c) A fc < T < A„. 

The analysis given above shows that the equation for the potential \i can be written in limiting cases as 

T 4 sin if - (Dr) a . b [dfi/dt - « 2 ] Q , 6 VQ a , b = (27) 
where the coefficients r a ^ and velocities v 2 b are equal to 

r „ - / ~v 2 v 2 { n/2 > T « Aa ' b f28) 
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FIG. 2: CM spectrum at intermediate temperatures (A& < T < A a ) for two cases: SI >> fi ai f, (a) and S7 << S~2 ra< b (b). There 
are two branches of the CM spectrum in the both limiting cases: the sound-like mode and the mode with a gap. 
The frequency and wave vector of CMs are normalized to S7* and fc* equal to: a) SI* = E\(T = 8Af,),fc* = Ei/vi and b) 
fl» = Ei(T — 2Ab), ft* = Ei/vi. With increasing temperature the two modes are getting closer to each other. 



where v 2 b = ^2(DA) a ^ b . The current density in each band is given by 



ja,b = (&p)a,bAa,bQa,b + <5a,fcE 



(29) 



with limiting values of conductivities 



, , i TP? - / exp(-A a)6 /T), T « A a , b 

(<rp)a, b - <?a,b \ nAa b/2 T : = a -' b \l, T » A Q , b 



(30) 



Qualitatively, dynamics of the CMs at any temperatures is described by Eqs. (|15ll8l27l29|) . One can exclude the 
varaibles Q Qjb and E and obtain the equations 



VaAa^a ~ V a B b {l b = -C Q (p, 



v b A b [i b - v h B a [i a - 1 



where A a = r a [-iQ + ( v a k) 2 {-itt + Q b )/M]; B a = ( i b k) 2 r b n a / M; M = ifi(if2-f2 + ); Q a = p a v a A a ; 
u a — Oajipa + &b) ■ Using the relation /i a — fit = {l/2)d^p/dt, we obtain the dispersion relation 



(31) 
(32) 

= Sl a + S7 b ; 



i£lv a v b (A a A b - B a B b ) = 2el(v a A a + v b A b - v b B a - v a B b ) (33) 

In the limiting cases T « A 0jb and A Qjb << T we obtain the formulas for the CM spectrum given above (see 
Ecis. (|21I25I26[) ). The "soft" , sound-like mode, which exists in the temperature range A a ^ b « T « T c , is of a special 
interest because another mode with a threshold frequency of the order of ea/h may be strongly damped if eo > Afc 
(this case seems to correspond to MgB 2 [ill]). 

The spectrum of this and other (plasma-like) mode may be obtained from Eq. (|33[) . Consider different limits of 
high and low frequencies Q with respect to the characteristic frequencies fl a , b - These frequencies in the considered 
temperature range (A b < T < A a ) are equal to: Q a = n(<T a /<T b )A a , il b = itA b /2T. 

cl) Assume that the frequency of oscillations is low 



n « si Qjb 



(34) 
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In the considered temperature range this condition implies that fl << A 2 /T << A;, << A Q . Then the functions 
A a fi and B ab acquire the form 

A a , b « (r B , 6 /in) [fi 2 - (k v a , b ) 2 ^} (35) 

B a , b « -Ka/iflJt^,,,) 2 ^ (36) 
Substituting these expressions for A a , b and -B a ,fc, we obtain the spectrum of oscillations 

n 2 = 2e\ + k 2 v\ (37) 

where e\ = eo[(j/ r ) _1 + (^h^fc) -1 ], u? = {v 2 Q b + wfn o )/0 +) r a = 1, r b = 7rA 2 /4T. We see that in this case the 
spectrum has a threshold frequency ~ (8 / ^e^T / (A b D b ) ■ The sound-like mode is strongly damped in this limit. 
c2) For large frequencies satisfying the condition 

Q » Sl a>b (38) 

we can write the functions A a j, and B ab in the form 

A a , b ^{r a:b /m)[n 2 -{kv a , b ) 2 }; B a , b «0 (39) 



The roots of Eq.j33J) are 

tf 1>2 = 2e 2 + \k 2 {i 2 a + i 2 ) ± {e\ + \k\v 2 a - i 2 ) + e 2 k 2 (i 2 a + i 2 ) 2e 2 k 2 {^ + -^-)}^ 2 (40) 
2 4 (rv) b {rv) a 

In the long-wave limit we obtain for fl 2 and fl 2 

fll = 2e\ + k 2 vl; ttl = k 2 vl (41) 

where v\ = {v 2 a (Dr) b + v\{vr) a j (v a r a + D b r b ), v 2 = {v 2 a (Dr) a + vl{Dr) b ) / (D a r a + D b r b ). Thus, we have a "hard" mode 
(f2i) and a "soft" mode (O2). However, the considered limit of high frequencies corresponds to a hypothetical case 
which is hardly realizable in experiments. We assumed that the frequency £1 is higher than £l ajb , but lower than A a ^ b . 
This means, in particular, that the inequality £l a < A& should be satisfied. This inequality can be presented in the 
form: cr a /a b << (Af,/7rA a ); that is, the conductivity of the a band with the higher energy gap A Q should be much 
less than the conductivity of the b band. 

In Fig. 2 we plot the dispersion curves fi(fe) for different temperatures determined by Eq. (l40|) . 



IV. OBSERVATION OF COLLECTIVE MODES IN A JOSEPHSON JUNCTION 

In this section we analyze an experimental method that allows one to observe the CMs and to determine the 
spectrum of these modes. The idea to identify the CMs in a superconductor has been suggested by Carlson and 
Goldman [1, [ijj and also used in Ref. • 

According to this idea one has to measure the I-V characteristics of a tunnel Josephson junction in the presence of 
a bias voltage Vb across the junction and of a weak external mag netic field H. As in Ref. [II], we consider a tunnel 
Josephson junction S\/ 1 / S2 consisting of two superconductors one of which is a two-band superconductor (S2) and 
another one (Si) is a single band superconductor. However, unlike Ref.(l6| we take into account a magnetic field, 
which allows one to determine the spatial dispersion of the CM spectrum and use the microscopic approach to derive 
necessary equations. As in Ref. [16j], we calculate the Josephson current jj taking into account terms of high order 
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in transparence |Tfc| 2 . As is well known, in the main approximation the Josephson current at a finite voltage Vb and 
in the presence of a magnetic field H has the form of a traveling wave (see, for example, [III [2?], Hl[ ) 

jj(x, t) = jj sm(n v t - K H x). (42) 

where fly = 2eVs/fi and Kh is given by Eq. (|62p . This current is injected into the superconductor S2 and modifies 
Eqs.([23| for the gauge-invariant potentials /i a ,b describing the conservation of charge of the Cooper pairs in each 
band. We show that, as in the Carlson and Goldman experiment [lH, a resonance occurs if the frequency f2 and the 
wave vector k coincide with Qy and Kh respectively. As was noted in Ref. 16], this resonance is analogous to the 
Fiske resonance in a Josephson tunnel junction when the velocity of the traveling wave (|42|) Vty /Kh coincides with 
velocity of the Josephson "plasma" waves. 

In order to generalize Eqs. ([23|) to the case of a Josephson junction, we assume that the thickness of the S2 electrode 
d is smaller than the London penetration depth Al as well as the quantity A: -1 and average Eq.(d]) over the thickness 



d taking into account the boundary conditions 24, 2 



(gdg/dz) a = (2R B a)- 1 [g s ^ a ] (43) 

where Rsa and a a are the interface resistance and conductivity for electrons in the a band in the normal state. The 
matrix Green's function g$ i n the one-band superconductor is assumed to describe a single band superconductor in 
equilibrium. The retarded (advanced) Green's functions in the S'-electrode have the standard form 

9s iA \z) = ih9s(e) + (if a cos Xs + ir x sin Xs)fs^)f {A) (44) 

with fg {A \e) = A 5 /£s (A) = {A s /e)gf (A \e), ^ {A) (e) = y/(e ± iO) 2 - A| ; X s is the phase in the S electrode. This 
averaging over the thickness d leads to the appearance in Eq. (fTOjl an additional term of the form iD a (2RB<j)~ 1 [gs,g a ]- 
The modified equation can be solved in the same way as it was done before. For example, the retarded function g^ 
is given again by Eq. (|12p with an additional term 



iEj a (g*-gXg*gX) (45) 

where the energy Ej a = D a j ' (2RB<j) a d is related to the interface transparency and to the Josephson coupling energy, 
d is the thickness of the two-band superconductor. This term results in a corresponding modification of Eq. (|27p . that 
acquires the form 

T 4 sin V - (vs)a,bdfi a ,b/dt + (PD/a) a , 6(Vj a ,h - (j a ,b/d) sin(xa,fc - Xs)) = (46) 

where the current j a .b at low temperatures is given by Eg ■ (| 1 7[) and the Josephson critical current j a ,b is equal to 
j a ,b — (A/ei?s) aj b, where the parameter A equals 

00 A A 

A a = 2nT -F== ~r==f ( 47 ) 



Cl 



At low temperatures A a = As ln(ciA a /A5) if As « A a and A a = A a ln(ciAs/A a ) if As >> A a , where 

= 2f™dxx\n(2x)/(x 2 + If/ 2 « 2.8. 

The phases Xa,b may be represented in the form 



Xa ,b = X±<P/2 (48) 

where x = (Xa + Xb)/2. 

We assume that there are no i^-solitons in the system (according to Ref. [J] special conditions are needed to create 
such solitons) so that in the main approximation one has: Xa = Xb] X ~ Xs = 2eVBt/h— Khx (see Appendix). 
Substituting these expressions into Eq. (|46|l . we obtain two equations 
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FIG. 3: Corrections to the I-V characteristics of a Josephson tunnel junction due to CMs in the two-band superconductor. 
The Josephson junction consists of a one-band and two-band superconductors. At voltages V satisying the equation Q v = 
f2§ + (KHVcm,) 2 a spike arises on the I-V characteristics. The position of the spike depends on an applied magnetic field H. 
The current is normalized to the value: jo = ^{jja — jjb){Ej a /i/ a — E 2 Jh jvh). 



-elsm(p-Uadn a /dt+(Pa)a^Qa = Ej a sin(fl v t - K H x), (49) 
eg sinip — v b d[i b jdt + (ra)i,VQi, = E 2 b sin(Sly£ — K H x), (50) 

where a a — 7r(crA) a /e and E 2 a — 2lS. a v a Ej a - Taking into account Eqs. (fT"5|) for Q a ,b, we can find the phase perturba- 
tion ip due to the "external forces" E 2 a b s'm(£lyt — Khx) 

<p = 2{Ej a /D a - EjjD^ImiN- 1 expi{fl v t - K H x)) (51) 
where N = fig — 57y(f7y + ij) + (KHV cm ) 2 . The phase perturbation <p leads to a change of the dc Josephson current 

= (jJa sm{n v t - K H x + S(x - Xs) + <p/2) + j Jb sm{n v t ~ K H x + S( X - Xs) ~ <p/2)) (52) 

8 X is the perturbation of the phase difference between the total phase x of the two-band superconductor and the phase 
Xs of the single band superconductor. The angle brackets mean the averaging in space and time. Then we expand 
the currents jj a with respect to S(x — Xs) and V 5 /2 and take into account Eq. (|51[) . The correction to the current jj 
due to the perturbation S( X — Xs) is omitted because it is not related to the CMs (this correction is real Fiske steps 
[H, HSHH) and corresponds to much larger voltages (see Appendix). The correction to the current jj due to <p is 

(Sjj) = -U/a " 3Jb)(Ej a /i?a - E 2 Ib /9 b )^ (53) 

where \N\ 2 = (ng - il 2 v + (K H v cm ) 2 ) 2 + ( 7 ft y ) 2 . 

We see that if the applied voltage satisfies the condition fly = ^lg+(KHV cm ) 2 , a spike arises on the I-V characteristics 
the position of which is shifted with varying magnetic field H . These spikes for different H are shown in Fig. 3. 



V. CONCLUSIONS 



In summary, using a simple model we have studied the CMs in two-band superconductors and shown that weakly 
damped CMs exist in these superconductors at all temperatures below T c in the frequency range: 7 < Q < A a ^ b , 
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where 7 is an effective charge imbalance relaxation rate. At low temperatures these modes have a spectrum similar 
to the spectrum of the Josephson "plasma modes" in tunnel Josephson junctions. However, the velocity of the CM 
in two-band superconductors by the order of magnitude is ~ %/ -DA, i.e., much smaller than the velocity of plasma 
modes in the Josephson junction. 

At high temperatures (T » A a ,&), the CM spectrum consists of two branches. One branch is analogous to the 
"plasma modes" in a Josephson junction but propagates much slower than the Swihart waves. Another branch is 
similar to the Carlson-Goldman mode in one-band ordinary superconductors, i. e., it has a sound-like spectrum. In 
the case of intermediate temperatures (A& < T < A a ) the gapless (sound-like) mode has a low damping if strict 
conditions are satisfied. 

The spectrum of the CMs considered above can be determined experimentally by measuring the I-V characteristics 
in a Josephson junction. Such a method has been applied by Carlson and Goldman [12] for measuring the CM 
spectrum in a conventional superconductor. They observed a small peak on the I-V curve the position of which 
depends on a weak applied magnetic field. In this case the wave vector k is proportional to H and the frequency of 
phase oscillations ft is related to the applied voltage: fl = (2e/H)VB- The measured temperature dependence of the 
spectrum of the CMs in two-band superconductors would allow one to elucidate the nature of superconductivity in 
such two-band superconductors as MgB%. 

Note that an evidence in favour of the existence of a Leggett-type collective mode in MgB2 was obtained in a recent 
work, where the point-contact spectrosopy was used [29| . 
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VII. APPENDIX 



Here we derive the relation between the phase difference 8 and the applied voltage V as well as the magnetic field 
H in the main approximation, i.e., in the absence of the Josephson coupling. In addition we present the derivation of 
Eq.(44). For simplicity we restrict ourselves with the case of low temperatures (T << A a> b). The results for higher 
temperatures are qualitatively the same. 

In the absence of the Josephson coupling the right-hand side of Eqs.(40) is equal to zero. Summing up these 
equations, we obtain the equation of the charge conservation 

eva^diiat/dt = {uD/a) atb Vj a ,b (54) 

where ja,b = ct aib (d x Xa,b — 2ttA x /&q)/2; $0 = hc/2e is the magnetic flux quantum. The coefficients a a ^ are related 
to London penetration depth Az,: a a ,b = c ^'o/(47r 2 A^ a b ). At low temperatures these coefficients are equal to a a> b = 
(aA) ai b/e. Summing up these equations, we obtain the equations of the charge conservation 

ed(v a fi a + D b ^b)/dt = V[a a Q a + a fc Q b ] (55) 

where Q a ,b = (dxXa,b ~ 2~kA x I$q)I2 is the condensate momentum. Differentiating this equation on time and taking 
into account Eq.(14) and the Poisson equation 

VE = Airp = 4Trve(D a (i a + ^bMb) (56) 

we obtain the equation 



v 2 d 2 {iy a p a + Vbp^/dt 2 = -k\ F {y a [i a + v h [ib) + V 2 [i5 a /z a + d b Hb] 



(57) 
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where v 2 = (7r/2)[(DpA) + (DPA) b ], v = v a + Vb, 8 a , b = {<TA) ayb /((aA) a + (crA) 6 ), fcf, F = 4ttcV The Thomas-Fermi 
screening length fc^, is of the order of interatomic spacing, i.e., is much shorter than charactetistic lengths of the 
problem: k TF >> v/A > v/ily ~ Kh- This means that the first term on the right-hand side in Eq. (T5T]) is much 
larger than other terms. Therefore we obtain 

VaV-a + v h [i h = (58) 

For the single band superconductor we have 

Ms = (59) 
Writing the potentials fi a ,b m the form \x a ^ = d(\ ± ip/2)/dt + eV, we obtain from Eqs. (|58H59|) 

d6/dt = 2eV -^{p a -D b )dip/dt (60) 

where 9 = \ ~ Xs'i Vb = V — Vs is the applied voltage. Eq. (|6"0|) generalizes the Josephson relation to the case of a 
junction with a two-band superconductor. 

If in the ground state there are no tp-solitons (ip = 0) and the external magnetic field H e is applied in the y-direction, 
one obtains from Eq. tfoT)]) 

e(x,t) = flvt-K H x (61) 
where fty = 2e Vb / h and the wave vector Kh is determined by the equation [27j 

K h = (27riJ e /$ ) [A tanh(d/A) + X s tanh(d 5 /As)] (62) 

where ds is the thickness of the single band superconductor, A = Xr, a + A^b- 

We turn now to finding solutions for Eqs. (|49H50|) . We assume that the phase perturbations are small (xa,b, <p « 1) 
and seek for a solution in the form 

(f(x, t) = Irrupu exp(£lyi — Khx) (63) 



Summing up Eqs. (|49H50j) . we obtain 

- d(P a fi a + D b n b )/dt + d[a a Q a + a b Q b ]/dx = Ej a + Ej b (64) 
The first term on the left is small due to the quasineutrality condition. Therefore we obtain from Eq. 

{a A) a 2 2 2 



(aA) b 

From the definition of Q a , b we have 



Q a + (Ej a + E 2 b )/{iK H v 2 D b ) (65) 



Qa ~ Qb = ^dep/dx (66) 



Finally we divide Eqs. (|49M50|) by v a , b and subtract from each other. We get 



+ ^)<P " \&<p/d? + 2v 2 a D a dQ a /dx = 2E 2 Ja ((,7! 



2 

Using Eq. (f6"6")l . we substitude Q a into Eq. (|6"7]) and come to Eo.([5*T|) for ip. 
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